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ABSTRACT 

Various computer enhancement and data extraction systems 
using LANDSAT data have been assessed and used to comple- 
ment a continuing geologic mapping program. Interactive 
digital classifioation techniques using both the parallelepiped 
and maximum -likelihood statistical approaches achieve very 
limited success in a> eas of highly dissected terrain. Com- 
puter enhanced imagery developed by color compositing 
stretched MSS ratio data has been constructed for a test 
site in northwestern Saudi Arabia. Initial results indicate 
that several igneous and sedimentary rock-types cAn be 
discriminated. 
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GEOLOGICAL MAPPING IN NORTHWESTERN SAUDI ARABIA 
USING LAND6AT-MULTISPECTRAL TECHNIQUES 


INTRODUCTION 

During the past 25 years, the Saudi Arabian Ministry of Petroleum and Mineral 
Resources, in cooperation «ith the U. S. Geological Survey, has actively en^ged 
in the geologic mapoing of the Precambiian areas (rf the Arabian Peninsula. Ini> 
tial field mapping efforts have contributed to the publication erf preliminary geo- 
logic maps at 1:500,000 and 1:2,000,000 (U.S. Geological Survey, 1963) scale 
for the entire peninsula, as well as a tectonic map at 1:4,000,000 (Brown, 1972). 

Present geolc^ic goals include mapping the entire shield at 1:100, 000 scale. The 
first few maps in this series have been published during the past several years, 
but a monumental task is still required to complete this ambitious project. 

Various approaches using LANDSAT data are currently being assessed and used 
to complement the continuing mapping program in this 600, 000 km^ area. The 
capability of discrimination of some of the lithologic units appears to be particu- 
larly promising. 

Figure 1 is a reduction made from a 1:1,000,000 mosaic constructed from 137 
LANDSAT frames. The entire Arabian shield is included in this regional cover- 
age, and the approximate margins of the Precambrian outcrop are delineated by 
the solid black lines. The basement rocks, in general, are seen as varying 
shades of dark gray. The very dark areas, on the other hand, are most com- 
monly the poorly reflective flood basalts of Tertiary or younger age. The flow 
identified on the figure provides an example of this poorly reflective rock class. 
The Precambrian outcrops extend approximately 2000 kilometers from the Gulf 
of Aqaba on the north, southward to the Gulf of Aden, and attain a maximum 
width of 700 kilometers, at approximately 23 degrees north latitude. 

Figure 2, an enlargement of the northera section of the mosaic, Includes the 
present area of interest. The outlined test sits area is located just east of the 
Red Sea coast in northwestern Saudi Arabia. 

Only limited and very generalized rock type distinctions are possible based only 
on photo interpretation of the various LANDSAT photographic formats. Several 
very generalized lithologic groupings can be distinguished on the basis of tonal, 
textural and distribution characteristics. These Include: (1) highly reflective 
alluvium and aeolian deposits, (2) poorly reflective, extensive basalt extrusions, 
(3) regionally homogeneous sedimentary strata and (4) coarse textured (at 
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llfure 1. Semi -font rolled mosaic of the Hed Sea and adjacent Precambrian 
basement. Original constructed from 137 LANDSAT scenes. 


2 



SINAt 




AEOLIAN 

'SEDIMENTS 


I.ANI)SA’I'-1 mosaic of northern IUhI Sea 
iind arliacent terrestrial areas. 




LANDSAT scale) massifs. Within the latter unit, sub-circular, highly reflective 
intrusions are common in the Red Sea area, and these, whore ma^^s are avail- 
able, are generally designated as being cf granitic composition. The high spec- 
tral reflectivity associated with these intruclons does not generally appear to be 
directly related to their discrete mineralogic composition, but rather is an in- 
direct indication of their topographic expression. Using the limited side-lap 
stereographic coverage available on LANDSAT imagery, light areas associated 
with these rocks are shown to be topographically low. This would suggest that 
the intrusicms are subject to more rapid mechanical weathering than their re- 
spective host rocks. The topographic depressions caused by such differential 
weathering appear now to be partially filled with highly reflective aeolian and/ 
or alluvial sediments. 

Identification of generalized rock-types such as those discussed above, however, 
are not generally adequate for most geologic mapping objectives. To improve 
on this discrimination capability, the use of LANDSAT digital data was evaluated 
using various computer enhancement techniques. 

The frame selected for the initial test site (1159-07272) is shown in Figure 3. 

This particular scene was selected because it includes two of the published 
1:100,000 geologic map areas. Thus larger scale geological mapping is avail- 
able to provide interim "ground truth" for about 20 percent of the scene outlined. 
Figure 4 includes those portions of the 1:500,000 published geologic maps of the 
northwestern Hljaz (Brown, Jackson, et al. , 1963) and northeastern Hijaz (Brown, 
Layne, et al. , 1963) which correspond to the test site. The slight shape varience 
of these identical geographic areas is caused by differences in image and carti- 
graphic projection. 

The geology included in the scene consists predominantly of complex Precambrian 
metamorphic assemblages intruded by plutonic rocks varying in composition from 
gabbro to granite. On the north, the Precambrian Is overlain by northeasterly- 
dipping sandstone of Cambrian (-€s) age. Quaternary /Tertiary flood basalts (QTb) 
make up the dark, less dissected areas in the west and southeast. Several small 
volcanic flows extend northward at the left side of the bottom margin. Because 
of the generally low spect.al reflectivity and poor tonal contrast, individual ig- 
neous and metamorphic units are difficult to distinguishln standard image formats. 

Figure 5 is a reproduction of the 1:100,000 Sahl al Matran (Hadley, 1973) map, 
which is entirely wlti.:.n the test site. The area included is shown in both Figures 
3 and 4. The additional detail provided by the larger scale is, of course, im- 
mediately recognized. In particular, it should be noted that the generalized 
granite units (gr) shown at 1:500,000 scale, have been divided and mapped as 
granite/granodiorite (ggd) on the south, and as two hornblende biotite quartz 
monzonltes intrusions (K>qm) on the north. These imits will subsequently be 
dit cussed in some detail. 
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APPHOACH 


In order to determine the applicability of computer classification and/or image 
enhancement for rook identification in this test ar ^a, the several programs avail- 
able at the Goddard Space Flight Center were as«.> ssed. Those included inter- 
active, multispectral classificaticm systems, employing both the maximum like- 
lihood and parallelepiped classification techniques and digital image enhancement 
using the SMIPS-VtCAR System*^ (Moik, 1973)* None of the systems have been 
exhaustively evaluated, but several preliminary conclusions can be projected. 


COMPUTER CLASSIFICATION 

Based on our preliminary studies ci lithology discrimination within the tost site, 
using interactive digital classification techniques, the more successful results 
were attained with the maximr.m likelihood approach incorporated into the Bendix 
M-DAS system. The most positive results attained with this system are shown 
on Figure 6, 

Training sites were selected for ten lithology classes, which fell into three gen- 
eral categories. These included: (1) Rocks mapped at 1:100,000 scale which did 
not appear distinctive in the color composite imagery. In this category, training 
sites were selected within the ITalaban group (hjm) and older greenstones (gd) 
mapped near the base of the Sahl «l Matran Quadrangle in order to determine if 
these units could be separated on the basis d their spectral differences. (2) Rocks 
somewhat distinctive on the composite imagery that correspond with field mapped 
units. These units included flood basalts (QTb) mapped on the Noi^hwest Hijaj 
Quadrangle and a granlte/granodiorlte (ggd) mapped at 1:100,000 on the Sahl al 
Matran Quadrangle. (3) Sediment units that appear to be distinctive on the com- 
posite imagery but have not been subdivided on the 1:100,000 geologic maps. 

The six training sites selected in this category were established in alluvium 
having different tonal characteristics on standard LANDSAT Imagery. In some 
cases, it appears that discrete alluvial units can be in part related to their spe- 
cific source rocks. 

Clearly, the units defined in red (basalt),** purple, yellow and dark blue (tonal 
alluvial units), correlate best with the geologic maps and the original MSS im- 
agery. The granite /granodiorite (ggd) though not uniquely distinctive, does ap- 
pear to be off-white, and its related alluvium depicted in dari: blue, defines the 
general outcrop area mapped on Figure 5, 


*501011 ln!eractivc Image Processing System, an interactive version of the Jet Propulsion Laboratory VICAR 
System. 

**AII color descriptions in this text arc based on original color imagery. Tonal changes arc generally visible on 
the ei)idvalent black-white figure format but these are considerably less definitive. 
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FlRure 3. NorthwesU>rn SiukU Arabia tost situ. Areas include' on Said al Matran 
(Hadley, 1973) and Wavban (Hadley, 1974) 1:100,000 ^colouic n aps are indicatetl 
by the northwest and southeast insets respectively. 
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Flinjro 4 Northwest Saudi Arabia test site neolonleal map. modified ;ifter Bnmn. 
.lacks.m. et al. (19 .3) and Brown. I.avne. et al. (19ti3». Orlulnal at 1 tf.OO. OOO 
scale. Areas of published 1:100,000 neolotrlc map coverage outlined. Note re- 
vised longitude sh'ft to west on larger scale maps. 
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Flj^rc 5, Geologic map of Sahl al Matran quadrangle, 
after Hadley (1973). Oriulnal at 1:100, ()0<) scale. 

On the other hand, no obvious discrimination can be made in the hinhlv dissected 
metamorphic and i^tneous terrains. It is important to note here that the more 
consistently correct classification Is associatiKl with the most homoj;eneous, 
near-horizontal beds formed by the flood basalts (red) and alluvium (purple, yel- 
low and dark blue). 

Fi({ures 7 and 8 are computer enlarKcments of portions of the l>asalt and granite/ 
Hrantxliorlte classes defined In Fin;ure 0. The small Imxes within each unit desi>{- 
nate the training sites selected for classification. In both fijjures, the nenerally 
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Figure 7. Computer enlargement of basalt (QTb) class using maximum-likelihood classi- 
fication. Angular box Is training site for this class. Note: Text based on original color 
imagery. Grey tonal chan,?es are less definitive. 
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Fi^jure 8. Computer enlargement of granite (ggd) class using maximum-likelihood classification. 
Rectangular box is training site for this class. Note: Text oased on original col jr imager>'. Grey 
tonal changes are less definitive. 





poor definition ci the class-unit margins is emphasized when compared with the 
equivalent mapped formations shown on Figures 4 and 5 respectively. 

Based on our preliminaty assessment of classification systems in this test site, 
it is concluded that the multispectral classification techniques appear to be se- 
verely limited for lithology mapping. One prime reason for this is the extreme 
difficulty encountered in defining adequate training sites in areas of highly dis- 
sected terrain. 


COMPUTER ENHANCEMENT 

Of the various rock discrimination approaches tested to date, a technique using 
color composites constructed from enhanced ratio data has proven to be most 
successful at this particular test site. An excellent summary of this method is 
provided by Rowan, et al. (1974). 

In simplified form, the procedure consists of three fundamental steps: (1) Selec- 
ted MSS bands are ratioed pixel by pixel, (2) the ratio data are stretched to en- 
hance contrast and converted into black-white positive film Imagery and (3) the 
stretched images of three ratio conibinations are combined into color composite 
imagery using various combinations of light filters and intensities. 

Ill the ratioing procedure, alternative ratio combinations are used to enhance 
specific rock classes. The 4/5, 5/6, and 6/7 combinations have provided the 
greatest single discrimination capability for both the Saudi Arabia test site and 
the Nevada test site of Rowan. In the Sahl al Matran area, however, other ratio 
combinations have successfully provided identification of additional rock classes. 
Optimum ratio combinations are often found to be scene dependent, and at present, 
no general rule exists for selecting ihe best co nbinations to provide the greatest 
discrimination capability over a wide range of climatic, topographic and/or light- 
ing conditions, Goetz and associates (Goetz, 1975, Goetz and Rowan, 1975) are 
conducting field studies using a portable field reflectance spectrometer to better 
define the optimum spectral and ratio combinations for specific applications. 

Both linear and nonlinear stretch processes can be alternatively used to optimally 
increase scene contrast over the dynamic range of the output product. Linear 
stretching Increases the scene contrast uniformly over the dynamic range while 
nonlinear stretching is used to increase local scene contrast in parts of the in- 
tensity range at the expense of other portions of the spectral range. In our test 
area, histogram equalization provided maximum intrascene contrast. Here 
again, as with the selection of ratio combinations, there are presently ho gen- 
eral rules that can be applied to all images to determine the optimum stretch 
parameters for a given set of cjnditlons. 
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We have been uiing a rotating drum film recorder (OPTRONICS) to conatruot 
black and white photo traniparenoiea from the digital imagea. The varicua 
oombinationa thua conatructed are then combined in a color -additive vie««r to 
rapidly determine the more aignificant oombinationa In regard to rock dikiorimi- 
nation capability. Once theae oombinationa are defined* they are prooeaaed 
through a CRT film recorder (DICOMED) to produce hard copy wit^ a better 
regiatrativin than that attained on the color-additive viewer. 

The reaultant compoaite providea the diatinot advantage of diaplaying a aingle 
type of material aa being aimil<'.r or the aame regardleaa of the local variation 
due to albedo and topography. Conversely* however, very diaaimilar materials 
easily separable on a standard composite image may beocnne inseparable on a 
ratio image where spectral reflectively slopes are similar. Thus* optimum 
rock-type discrimination must be accomplished in several combinations of dif- 
ferent ratio oombinatians. 

All digital enhancements described in this study were performed using an inter- 
active image processing system on IBM/360-75 or -91 computers. 

Figure 9 is a DICOMED ratio-composite where the MSS band 4/5 ratio im- 
age is projected with a blue niter* 5/6 image with a green filter* and 6/7 
image with a red filter. Note in this presentation that the monzonites (hbqm) 
mapped on the Sahl al Matran Quadrangle are a dark green while the granite/ 
granodiorite (ggd) and alluvium around the granites are ocher in hue. Outside 
the Matran map area, at the right margin, a dark green triangular body has also 
been mapped as a homblende-biotite quartz monzonite on the 1:100*000 Wayban 
Quadrangle (Hadley* 1974), In addition the ocher-brown body Immediately to 
the west of the monzonite is mapped as rlebeckite granite (rg). There is no 
1:100, 000 coverage of the ocher-colored area to the north, but it is also shown 
as granite (gr) on the 1:500,000 Northwestern Hijaz Quadrangle, 

The brilliant blue* unit on the composite, of course* corresponds to the basalt 
flow (QTb) classified on Figure 6. It can be seen very clearly here that the mar- 
gins of the basalt are much more distinctly defined than they were using a com- 
puter classification system (Fig. 7), In addition to these flood basalts, several 
smaller flows are seen near the eastern margin of the image and several dis- 
crete cinder cone fields are visible across the wadi from the large flow. This 
enhancement has been used to aid in distinguishing small cinder cones of basaltic 
composition from ultramafic rocks in the preliminary 1:100,000 map compilntion 
for the area south of the Sahl al Matran Quadrangle. 


•Blue reproduce* as light grey in the black-white format and therefore it is difRcult to distinguish basalt 
from the underlying schist in Figure 9. Please refer to Figures 1 1 and 12 for better definition of these flows. 
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Figure 9. DICOMED ratio-composite of a part of the northwest Saudi Arabia test site. Constructed 
from MSS 4/5 ratio imagery projected through a blue filter, MSS 5 G imagery through a green filter 
and MSS 6/7 imagery through a red filter. Note: Text based on original color imagery. Grey tonal 
changes are less definitive. Refer to Figures 11 and 12 for better definition of basalt. 


The tan-colored outcrops in this enhancement are mostly complex metamorphic 
assemblages. Our ability to distinguish between metamorphic units mapped in 
the field has been poor to date. 

Figures 10 and 11 were photographed from a color additive viewer screen. Some 
resolution degradation caused the poor registration of the different ratio images 
results in a slight blurring. The points of reference, however, can still be 
readily distinguished. 

The Figure 10 ini!^;e was made from a combination using MSS band t-atio 4/7 
imagery with blue filter, 7/4 imagery with a green filter, and 6/5 imagery with 
a red filter. In thi:^ renditicm, the Siq sandstone (^es) in the upper right, is a 
bright yellow. This is not unique to the scene, however, for granites (gr) to 
the south are of the same hue. It is thus not possible to distinguish between the 
two different lithologic units on the basis of this particular ratio enhancement 
alone. When the MSS 6/5 band ratio image is replaced by 7/6 image using a red 
filter, however, tonal changes occur as shown in Figure 11. In this second ratio 
composite, the Siq sandstone remains a bright yellow, but the granites that were 
yellow in Figure 10 are now a pale yellow-green hue. Conversly, the monzonlte 
(hbqm), that is pale yellow-green In Figure 10, is now yellow on Figure 11 with 
the new ratio combination. The use of a pair of ratio combinations, thus, seems 
adequate in this case to define the S^q sandstone within this particular scene; at 
the same time, the inherent spectral characteristic differences between granites 
and monzonltes are reemphasized. 

Figure 12 is a DICOMED ratio-enhancement that overlaps the area shown in 
Figure 9, on the south. This color ratio-composite was made using MSS 5/4 
imagery with a blue filter, 6/5 Imagery with a green filter, and 7/6 Imagery 
with red filter, or just the inverted ratio combinations used in the construction 
of Figure 9, The several volcanic flows that extend northward from the left 
side of the bottom margin of the test site (Fig. 3) are shown in two colors in this 
ratio-enhancemen';. Based on field confirmation, the geologically recent flows 
correlate with the areas shown in ocher while the flows of historic age are shown 
as black (the most recent eruptions in the area were recorded in 1280 A. D. ). 

The flood basalts (QTb) that were blue in Figure 9 are a bright yellow In Figure 
12. There is, thus, a strong suggestion that the progression of color tonal dif- 
ferences from black to ocher to yellow among the basalts can be at least crudely 
related to age. Fundamental differences in chemical composition of tne respec- 
tive flows could provide an explanation for the color differences. Field data, 
however, suggest that the compositi onal differences are not sufficient to signifi- 
cantly alter the spectral response between the different volcanic units. Changes 
in reflectance due to a varying degree of weathering, on the other hand, better 
explain the observed difference of spectral response, and of course, this would 
also relate to age differences. 
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Figure 10. Ratio-composite enhancement photographe<i from a color additive viewer screen. 
Composite is constructed from the MSS 4 7 ratio image projected through a blue filter, the 
MSS 7 4 image through a green filter and MSS the fi 5 image through a r ’ filter. Note: Tes 
based on original color imagerv. Grev tonal changes are less definitive. 



Fifjurc 11. Ratio-composite enhancement photORrapheti from a color-additive viewer screen. 
Composite is constructe<i from MSS 4/7 ratio imagery projecteti through a blue filter, MSS 
7 '4 imagery' through a green filter ami MSS 7 'f. imagery through a reti filter. Note: Text 
b.'ise<l on original color imagery. Grey tonal changes ar-? less definitive. 




In addition to tho mapped extrusivoa dlsoussod above, thero appears to be a dis- 
tinct lithologic class defined in pale yellow in the southeast comer of Figure 12, 
The entire lithologic suite in this quadrant has been mapped as schist (so) at 
small scale (see Fig. 4), Tho enhanced imagery was compared with a 1:100, 000 
orthophoto mosaic the area which had been constructed from reduced 1:30, 000 
aircraft photography, and it was found that the pale-yellow spectral unit defined 
on the ratio-enhanoed image correlated extremely well with a clearly defined, 
dark, "botryoidal"-textured unit on the mosaic. It is anticipated that this unit 
v'lU bo investigated in tho field to determine if tho pale-yellow unit does indeed 
define a specific metavolcanic sequence. 

Most of the western one-fourth of Figure 12 is within the area to be included in 
the 1:100, 000 Dhaylan Quadrangle. The lithology change across tho well-defined 
north-south trending fault (A-A*) Is especially well displayed on tho south where 
an abrupt color change occurs from purple on the west to a dark-brownish-green 
on the ^st. The southern part of the pale green sinuous body (gr) is designated 
as granite on the preliminary geologic map. On the same map, however, this 
granite is shown to terminate at a secomlary fault extending between B and B* on 
Figure 12. The ratio composite data, however, strongly suggest that the granite 
continues in a northwesterly direction as a discrete unit, and terminates at the 
major north-south fault. This area must be further Investigated in the field to 
determine the true extent of the granite body. 

To complement the spectral-lithologic relationships discussed above, special 
attention was placed on the areas immediately adjacent to the three ancient mine 
workings described by Hadley (1972). Each area was studied in detail using a 
variety of ratio/filter and stretch combinations in order to identify any hydro- 
thermal alteration that might be associated with the mineralization. No anomalies 
were recognized, but this could very well be due either to the absence of altera- 
tion products or the small areal extent of mineralized areas. Plans are presently 
being formulated to study several well-developed hydrothermal alteration zones 
known to exist in association with the massive Wadi Wnssat sulfide deposits near 
the Yemen border. 


CONCLUSIONS 

The ability to identify lithologic classes using LANDSAT multls|>ectral techniques 
appears to be possible, at least to a limited degree. The various classification 
techniques that have been so successful in a variety of agriculture, land use and 
environmental studies appear to be severely limited in rock-type discrimination, 
especially in areas of highly dissected terrain. This is not necessarily the re- 
sult of the inadequacy the various statistical approaches used in classification, 
but rather, the extreme difficulty encountered in defining adequate training sites 
in areas of highly dissected terrain. 
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A technique of image enhanooment using stretched ratio-oomposite data uses 
slope difference as the basis for discrimination rather than absolute reflectance 
values. By using this approach, spectral differences caused by the effects of 
local variation of terrain and albedo are minimissed. 

Preliminary studies, using images constructed using this technique, indicate 
that several rock-types can be distinguished. Those include basaltr, mcmzonites, 
and certain granites. Efforts to distinguish between different metamorphio lithol- 
ogies have been unsuccessful. 

Computer enhanced LAND6AT imagery continues to be used in selected tost sites 
to complement traditional mapping techniques in Saudi Arabia, The initial prom- 
ising work is being refined and expanded to form a base from which we can ex- 
pand and refine our present capabilities. 
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